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Using a Maximum Entropy Production Principle (MEPP), we derive a new type of relaxation 
equations for two-dimensional turbulent flows in the case where a prior vorticity distribution is 
prescribed instead of the Casimir constraints [Ellis, Haven, Turkington, Nonlin., 15, 239 (2002)]. 
The particular case of a Gaussian prior is specifically treated in connection to minimum enstrophy 
states and Fofonoff flows. These relaxation equations are compared with other relaxation equations 
proposed by Robert & Sommeria [Phys. Rev. Lett. 69, 2776 (1992)] and Chavanis [Physica D, 
237, 1998 (2008)]. They can serve as numerical algorithms to compute maximum entropy states 
and minimum enstrophy states with appropriate constraints. We perform numerical simulations of 
these relaxation equations in order to illustrate geometry induced phase transitions in geophysical 
flows. 
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I. INTRODUCTION 

Two-dimensional incompressible and inviscid flows are 
described by the 2D Euler equations 



u- Vu; = 0, 



u = — z x V^, (1) 



where ujz = V x u is the vorticity, ip the stream func- 
tion and u the velocity field (z is a unit vector nor- 
mal to the flow). The 2D Euler equations are known to 
develop a complicated mixing process which ultimately 
leads to the emergence of large-scale coherent structures 
like jets and vortices |l|-|4( . The jovian atmosphere shows 
a wide diversity of coherent structures [5|-[9j: Jupiter's 
great red spot, white ovals, brown barges, zonal jets,... 
Similarly, in the earth atmosphere and in the oceans, 
there exists large-scale vortices such as modons (pairs of 
cyclones/anticyclones) or currents like the Gulf Stream or 
the Kuroshio Current. One question of fundamental in- 
terest is to understand and predict the structure and the 
stability of these quasi stationary states (QSSs). This 
can be done by using elements of statistical mechanics 
adapted to the 2D Euler equation. 

In recent years, two statistical theories of 2D turbulent 
flows have been proposed by Miller-Robert- Sommeria 
(MRS) and Ellis-Haven-Turkington (EHT). These the- 
ories mainly differ in the type of constraints to be con- 
sidered. Miller [Toj| and Robert & Sommeria [ll| assume 
a purely conservative evolution (no forcing and no dissi- 



pation) and take into account all the constraints of the 
2D Euler equations. The equilibrium state is obtained 
by maximizing a mixing entropy S[p] while conserving 
energy, circulation and all the Casimirs. On the other 
hand, Ellis, Haven & Turkington [l2| argue that, in real 
flows undergoing a permanent forcing and dissipation, 
some constraints are destroyed. They propose a prag- 
matic approach where only the robust constraints (energy 
and circulation) are taken into account while the fragile 
constraints (Casimirs) are treated canonically. This is 
equivalent to introducing a prior vorticity distribution 
x(cr). This prior vorticity distribution is assumed to be 
determined by the properties of forcing and dissipation. 
The equilibrium state is then obtained by maximizing a 
relative entropy S x [p] (depending on the prior) while con- 
serving only energy and circulation (robust constraints). 

Some relaxation equations towards the statistical equi- 
librium state have been proposed in each case. Consider- 
ing the MRS approach, Robert & Sommeria [l3| obtained 
a relaxation equation from a Maximum Entropy Produc- 
tion Principle (MEPP) by maximizing the production of 
entropy S[p] at fixed energy, circulation and Casimirs. 
On the other hand, considering the EHT approach, Cha- 
vanis 

[313 

showed that the prior determines a general- 
ized entropy S\uJ] and that the mean flow is a maximum 
of generalized entropy at fixed energy and circulation. 
He then obtained a relaxation equation from a MEPP by 
maximizing the production of generalized entropy S[uJ] at 
fixed energy and circulation. Interestingly, the resulting 
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equation has the form of a nonlinear mean field Fokker- 
Planck (NFP) equation that appears in other domains of 
physics Hzllift 

In this paper, we introduce a new class of relaxation 
equations associated with the EHT approach by maxi- 
mizing the production of relative entropy S x [p] at fixed 
energy and circulation in order to obtain the evolution 
of the full distribution of vorticity levels. Interestingly, 
this leads to a new class of relaxation equations that does 
not appear to have been studied so far. We derive the 
corresponding hierarchy of moment equations and show 
that it is closed in the case of a Gaussian prior leading 
to a minimum enstrophy state. These relaxation equa- 
tions can serve as numerical algorithms to compute max- 
imum entropy states with appropriate constraints. In the 
present paper, we develop the theory and discuss in de- 
tail the link between the MRS and the EHT approaches. 
We also give a numerical illustration of our relaxation 
equations in relation to minimum enstrophy states and 
Fofonoff flows. 

The paper is organized as follows. In Sees. HH IIIII and 
II V Al we provide a short review and comparison of the dif- 
ferent relaxation equations introduced in the context of 
2D turbulence in relation to the MRS and EHT statistical 
theories. This review can be useful to people interested 
in this topic and should clarify the subtle connections 
between the different equations. We also size this oppor- 
tunity to improve the discussion and the presentation of 
these relaxation equations and explicitly treat particular 
cases. In Sec. IIVB[ we introduce a new class of re- 
laxation equations adapted to the EHT approach based 
on the specification of a prior vorticity distribution. In 
Sec. [VI we consider the case of a Gaussian prior where 
these equations can be simplified. In Sec. IVI[ we use 
these relaxation equations to illustrate phase transitions 
in geophysical flows, in particular the transitions between 
monopoles and dipoles when the domain becomes suffi- 
ciently stretched [ll-Elj. 

Note: in this paper, we shall mainly follow the pre- 
sentation of Ellis et al. [l2j] and Chavanis [lil-HiSj who 
first introduced the maximization problems ([36]) and ([52]) 
based on the notion of priors. However, in Sec. IIII CI 
we note that these maximixation problems also provide 
sufficient conditions of MRS thermodynamical stability 
(this is t he p resentation adopted by Bouchet [22] and 
Chavanis [23j). Therefore, the study of the maximization 
problems ([36]) and ([52]) , and the corresponding relaxation 
equations, is interesting in these two perspectives. 



II. STATISTICAL MECHANICS OF VIOLENT 
RELAXATION 

A. The Miller- Robert- Sommeria theory 

Starting from a generically unsteady or unstable ini- 
tial condition, the 2D Euler equation develops an intri- 
cate filamentation leading, on the coarse-grained scale, 



to a quasi stationary state (QSS). The problem is to 
predict the structure of this QSS as a function of the 
initial condition. A statistical theory of the 2D Euler 
equation has been proposed by Miller [To| and Robert & 
Sommeria pjj (see also Kuz'min (24|) by extending the 
approach of Onsager based on point vortices [25l-[27( or 
the approach of Kraichnan based on the truncated Euler 
equations [H, HI- The Miller-Robert-Sommeria (MRS) 
theory is analogous to the theory of violent relaxation 
developed by Lynden-Bell [30] in the case of collisionless 
stellar systems described by the Vlasov-Poisson system 
(see Chavanis 0, [3l| for a description of the numerous 
analogies between 2D vortices and stellar systems). The 
key idea is to replace the deterministic description of the 
flow u)(r,t) by a probabilistic description where p(r, cr, t) 
gives the density probability of finding the vorticity level 
uj = a in r at time t. The observed (coarse-grained) 
vorticity field is then expressed as cJ(r, t) = J pcrdcr. The 
MRS approach is well-suited to isolated systems (no forc- 
ing and no dissipation) where all the inviscid invariants of 
the 2D Euler equation are conserved. The MRS statisti- 
cal equilibrium state is obtained by maximizing a mixing 
entropy 



S[p] = — J p(r, cr) In p(r, cr) drdcr, 



(2) 



while respecting the normalization condition J pda = 1 
and conserving all the inviscid invariants which are the 
energy 



E = 



1 



uo^dv 



the circulation 



r = 



^pcr drda, 



pa drdcr, 



(3) 



(4) 



and the Casimirs invariants If = J f(u)dr. This includes 
in particular the conservation of all the microscopic mo- 
ments of the vorticity 



_ / — , 

■ n>l 



/""*-/ 



pa n drdcr. 



(5) 



The conservation of the Casimirs is also equivalent to 
the conservation of the fine-grained vorticity distribution 
7(a) = J p(r, a) dr, i.e. the total area j(a)da occupied 
by each vorticity level cr. 

We shall distinguish between robust and fragile con- 
straints. This distinction will be very important in the 
following. The circulation and the energy are called ro- 
bust constraints because they c an be expressed i n term s 
of the coarse-grained vorticity: T[lj] ~ T[uj] and E[uo) ~ 
E\lu] (the energy of the fluctuations can be neglected 
[IlJuJj)- By contrast, the higher moments of the vortic- 
ity are called fragile constraints because, when calculated 
with the coarse-grained vorti city cJ, they are not con- 
served since cu™ ^ uJ n . Thus r n> i[o;] ^ r n> i[cJ]. They 
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must be expressed therefore in terms of the "microscopic" 
vorticity distribution p(r, a) as in Eq. (|5]) where we have 
introduced the local moments uo n = J per 71 da of the vor- 
ticity distribution. We must therefore distinguish the mi- 
croscopic moments of the vorticity r{'| 1 [p] = J uj n dr — 
J pa n drda (conserved) from the macroscopic moments 
of the vorticity r^'-JcJ] = J uJ n dr (non-conserved). 

In the MRS approach which takes into account all the 
constraints of the 2D Euler equation, the statistical equi- 
librium state is determined by the maximization problem 



max {S[p\ | Tp] = T, Ep] = E, 

p 

r£|i[p] = r£»i, j P da = i}. 



(6) 



The critical points of mixing entropy at fixed E, T, 

[53| and normalization are obtained from the variational 

principle 

5S-f35E-a5r-J2 a nSF f n 9 -- J ((r)5pdadr = 0, (7) 



n>l 



where /?, a, a n> \ and £(r) are appropriate Lagrange mul- 
tipliers. This leads to the Gibbs state 



P0,cr) 



1 



Z(r) 



(8) 



where Z = J x(<j)e _ ^^ +a ^ )cr (i(j is a normalization 
factor and we have defined the function x( cr ) = 
exp(— Y^ n >i a n° n ) which encapsulates the Lagrange 
multipliers associated with the fragile constraints. The 
coarse-grained vorticity is then given by 



uj = 



1 dlnZ 



Jx(cr)e-^^ a da ' (3 # 
The function F is explicitly given by 
F($) = -(Inx)'(*), 



= F(M + a). (9) 



(10) 



where we have defined x($) = / x( cr ) e_cr< ^^ cr - O n the 
other hand, differentiating Eq. ^ with respect to -0, it 
is easy to show that the local centered variance of the 
vorticity distribution 



0J2 



■uj 2 



p(cr — uj) 2 dr, 



is given by 



0J2 



1 1 d 2 lnZ 



(11) 



(12) 



As noted in [9], this relation bears some formal simi- 
larities with the fluctuation-dissipation theorem (FDT). 
Since uj = cJ(V0? the statistical theory predicts that 
the coarse-grained vorticity cJ(r) is a stationary solu- 
tion of the 2D Euler equation. On the other hand, since 



uJ f (ip) = —f3uj2(^) with UJ2 > 0, the uj — ip relationship is 
a monotonic function that is increasing at negative tem- 
peratures f3 < and decreasing at positive temperatures 
f3 > 0. Therefore, the statistical theory predicts that the 
QSS is characterized by a monotonic uj(ip) relationship. 
This uj — ^ relationship can take different shapes depend- 
ing on the initial condition. Substituting Eq. (J9j) in the 
Poisson equation (pQ), the equilibrium state is obtained 
by solving the differential equation 



l_d_ 

~]3~ckp 



In Z, 



(13) 



with ijj = on the boundary of the domain and relating 
the Lagrange multipliers to the constraints. Then, we 
have to make sure that the distribution (j8j) is a (local) 
maximum of entropy, not a minimum or a saddle point. A 
critical point of constrained entropy is a (local) maximum 
iff M : 



S 2 J 



\J5p I SE 



(5pf 
P 

ST = ST 



drda 



7} J SuSiIj 



dr < 0, 
da = 0, (14) 



i.e. for all perturbations Sp that conserve the constraints 
(circulation, energy, Casimirs, normalization) at first or- 
der. Finally, if several (local) entropy maxima remain for 
the same values of the constraints, we can compare their 
entropies to determine which one is the global entropy 
maximum and which one is a relative entropy maximum. 
We stress, however, that local entropy maxima can be 
long-lived, hence fully relevant, for systems with long- 
range interactions [Hi). 



B. Relaxation equations 

Robert & Sommeria [13| have introduced a relaxation 
equation solving the optimization problem (J6j) by maxi- 
mizing the rate of entropy production S at fixed circula- 
tion, energy and Casimir constraints (and other physical 
constraints putting a bound on the diffusion currents). 
This Maximum Entropy Production Principle (MEPP) 
leads to the following relaxation equation 

^ + u • Vp = V • [D (Vp + p(t)(<r - 5J)pW)] , (15) 



where 



f3(t) 



J DVZJ-Vipdr 
J £>w 2 (W>) 2 dr' 



(16) 



is a Lagrange multiplier (inverse temperature) enforcing 
the energy constraint E = at any time and D(r,t) is 
a diffusion coefficient. The diffusion coefficient is not de- 
termined by the MEPP but it must be positive to have 
an increase of entropy (see below). The boundary condi- 
tions are J n = where J = —D(Vp J r f3(t)(a — cJ)pV?/0 



is the current of level a and n is a unit vector normal to 
the boundary. Easy calculations lead to the H-theorem 



S= I ^ (Vp + P(t)p(a - uj)V^f drda > 0. (17) 



Equation (fT5j) with the constraint ([15]) has the following 
properties: (i) T, E, and j pda = 1 are conserved, 

(ii) S > 0. (hi) S = ^ p(r,cr) is the Gibbs state © 
<^> c^p = 0. (iv) p(r, a) is a steady state of Eqs. ([T5]) - 
(p~6|) iff it is a critical point of S at fixed E, T, r{'| 1 and 
normalization, (v) A steady state of Eqs. (fT5]) - ([T6|) is lin- 
early stable iff it is a (local) maximum of S at fixed E, T, 
r{'| 1 and normalization. By Lyapunov's direct method, 
we know that if S is bounded from above, Eqs. (fT5]) - (fT6|) 
will relax towards a (local) maximum of S at fixed E, T, 
and normalization (if several local maxima exist, 
the choice of the maximum will depend on a notion of 
basin of attraction). Therefore, a stable steady state of 
Eqs. (fT5]) - (fT6|) solves the maximization problem (j6j). By 
construction, it is the MRS statistical equilibrium state 
(most mixed state) corresponding to a given initial condi- 
tion. As a result, the relaxation equations ([T5|) - (fT6]) can 
serve as a numerical algorithm to solve the maximization 
problem (|6j) for a given value of the constraints specified 
by the initial condition. These relaxation equations have 
been studied theoretically and numerically in (33l436|. 



C. Moment equations 



Remark 1: the relaxation equations (|T5]) - ([T6|) do not 
respect the invariance properties of the 2D Euler equa- 
tion (invariance by translation or rotation) and this may 
be a serious problem to describe the evolution of the flow 
into several isolated vortices (if we use these equations 
as a parametrization of 2D turbulence). A solution to 
this problem has been proposed by Chavanis & Somme- 
ria [37j by reformulating the MEPP under a local form, 
introducing currents of energy, angular momentum and 
impulse. However, the resulting relaxation equations are 
more complicated and have not been numerically solved 
for the moment. 

Remark 2: in order to take into account incomplete 
relaxation [H, [39| , Robert & Rosier [33| and Chavanis 
et al. [3l| have proposed to use a diffusion coefficient 
depending on the local fluctuations of vorticity. This can 
freeze the system in a "maximum entropy bubble" which 
is a restricted maximum entropy state [38| . This diffusion 
coefficient can also be justified from a quasilinear theory 
of the 2D Euler equation [40]. It can be written in the 
form Ella!: 



D(r,t) 



2 1/2 



(18) 



where K is a constant of order unity, e the scale of unre- 
solved fluctuations and 002 the local centered variance of 
the vorticity. 



From Eq. ([T5]) , we can derive a hierarchy of equations 
for the local moments of the vorticity uo n = J per 71 da. 
The equation for the moment of order n is 



duj n 



V <D 



Vuu n + f3(t)(uj n+l - uJu n )V^ 



(19) 



This hierarchy of equations is not closed since the equa- 
tion for the moment of order n involves the moment of 
order n + 1. Robert & Rosier [33| have proposed to 
close the hierarchy of equations by assuming that the 
density distribution p(r, a, t) maximizes the entropy (j2j) 
with the constraints of the known first n moments and 
the normalization. This gives a density of the form 
p ~ exp(— Y^k=i a k°~ k ), where the Lagrange multipliers 
at can be calcul ated fr om the constraints of the known 
moments. Then, o; n+1 can be obtained from this distri- 
bution and expressed in terms of cJ,...,^. 

Kazantsev et al [4l| have considered in detail the case 
n = 2. If we maximize the entropy ([2]) at fixed uJ(r, t) = 
J pa da and o; 2 (r,t) = J pa 2 da, we obtain a Gaussian 
distribution of the form 



1 



(20) 



Therefore, the vorticity distribution is locally Gaussian 
with mean value cJ(r, t) and centered variance c^2(r, t). 
From this distribution, we compute u 3 = u 3 + 3ujuj2 
which closes the hierarchy at the order n = 2. Then, 
the equations for the mean vorticity and the centered 
variance can be conveniently written 



— + u • VcJ = V 

at 



D 



VcJ + p{t)uj 2 V^ 



(21) 



duj2 
~dt 



U • Vco>2 



V • (DVuj 2 ) 

2DVuj • (VcJ + /3(£)cj2 V?/0 • (22) 



These equations conserve by construction the energy, the 
circulation and the microscopic enstrophy T^' 9 ' = J oo 2 dr 
but not the higher moments. On the other hand, we can 
prove an i^-theorem for the entropy. Using the Gaussian 
distribution (|20]h the mixing entropy (|2]) can be written 



\n0J2 dr, 



(23) 



up to some unimportant additive constants. Then, using 
Eqs. ([2T]) and ([22]) . it is easy to establish the H-theorem 



dr+ / — (VcJ + /3(t)cj 2 V^) 2 ^r > 0. 

UJ2 



(24) 
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At equilibrium, S = 0, we obtain 

VcJ + (3uj 2 0, 

Vo; 2 0. 



(25) 



(26) 



The second equation shows that the centered variance 
Co>2 (r) = Q2 is uniform and the first equation can then be 
integrated into 



uj = -(] 2 (^ + a), 



(27) 



where a is a constant. The statistical equilibrium state 
presents a mean flow characterized by a linear uJ — tp 
relationship and Gaussian fluctuations around it. It is 
a maximum of entropy at fixed energy E, circulation T 
and microscopic enstrophy T^' 9 ' . It is also a minimum of 
macroscopic enstrophy T^' 9 ' at fixed energy E and circu- 
lation T (see Appendix [XJ) . 

Remark: We note that, with this formalism, it is tech- 
nically difficult to go beyond the Gaussian closure ap- 
proximation. The approach developed in Sec. IIV Al may 
provide an alternative strategy to describe more complex 
situations where the vorticity distribution is not Gaus- 
sian. 



D. The equation for the velocity field 

Following Chavanis & Sommeria 0|, we can derive a 
relaxation equation for the velocity field. The equation 
for the coarse-grained vorticity field is 



duj 
~dt 



u • VcJ = V • < D 



(28) 



where we recall that D(r 1 t) can depend on position and 
time. For a 2D field, we have the identity V x (z x a) = 
(V-a)z. Therefore, we can rewrite the foregoing equation 

as 



^ + u • VcJ 



(29) 



Since V x u = cjz, the corresponding equation for the 
velocity field is 

+ (u • V)u = --Vp + z x D (Vw + P(t)w 2 Vi/>) , 



(30) 



where p is the pressure and p the density. Now, using 
u = — z x and the identity 



Au = V(V • u) - V x (V x u) 

= -V x (cJz) = z x VcJ, 



valid for a 2D incompressible flow, we finally obtain 

^ + (u • V)u = --Vp + D (Au - (3(t)u 2 u) . (32) 
at p 

We see that the drift term in the equation for the vorticity 
(f28|) takes the form of a friction in the equation for the 
velocity ([32]) . Furthermore, the drift coefficient or the 
friction coefficient is given by an Einstein-like formula 
£ = D/3 involving the diffusion coefficient and the inverse 
temperature. At equilibrium, we get 



Au = /Jtt^u, 



(33) 



which can be directly derived from the Gibbs state (j8j) 
using Au = z x VcJ, VcJ aJ'(^)V^ and Eq. (fT2j) . 



III. STATISTICAL MECHANICS WITH A 
PRIOR VORTICITY DISTRIBUTION 

A. The Ellis-Haven- Turkington approach 

In the MRS theory, it is assumed that the flow is rig- 
orously described by the 2D Euler equation so that all 
the Casimirs are conserved. However, in many geophys- 
ical situations, the flows are forced and dissipated at 
small scales (e.g., due to convection in the jovian atmo- 
sphere) so that the conservation of the fragile constraints 
(Casimirs) is destroyed [?J. Ellis, Haven and Turkington 
[12| have proposed to treat these situations by keeping 
only the robust constraints E and T and replacing the 
conservation of the fragile constraints by the speci- 
fication of a prior vorticity distribution x( cr )- As noted by 
Chavanis [15[ , this amounts to making a Legendre trans- 
form of the MRS entropy (j2j) with respect to the fragile 
constraints (see Eq. ([34]) ). The EHT approach corre- 
sponds therefore to a grand microcanonical version [23| 
of the MRS theory in which the chemical potentials a n >i 
associated with the microscopic constraints are given [55j. 

We introduce the grand entropy [l5| : 



s x [p] = %]-E a » r « ! 



n>l 



Explicitly, we have 



S x [p] 



J p(r,a) In 



P(r,c) 



drdcr, 



(34) 



(35) 



(31) 



where x( cr ) = exp(— J2 n >i a n^ n )- In the present con- 
text, this function is given and is called the prior vorticity 
distribution. It is determined by the properties of forcing 
and dissipation for the situation considered. On the other 
hand, S x is called the relative entropy [7]. The EHT sta- 
tistical equilibrium state is obtained by maximizing the 
relative (or grand) entropy (|35|) while respecting the nor- 
malization condition J pdcr = 1 and conserving only the 
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robust constraints r = J uodr and E = - J ujipdr. There- 
fore, we have to solve the maximization problem [l2j: 



max {S x [p] | r[ZJ] = T, £[cJ] = E, J pda = 1}. 



(36) 



The critical points of grand entropy S x at fixed cir- 
culation, energy and normalization (canceling the first 
variations) are given by the variational principle 

SS X - f35E -aST- J ((r)5pdadr = 0. (37) 

This leads to the Gibbs state flSJ. Therefore, the critical 
points of the variational principles (|6j) and ([36]) coincide. 
On the other hand, a critical point of constrained grand 
entropy is a (local) maximum iff (23| : 



2 T - 



S 2 J 



{5pf 



drda 



Suj5ip dr < 0, 



V5p | 5E = ST = J Spda = 0, 



(38) 



i.e., for all perturbations dp that conserve circulation, en- 
ergy and normalization at first order. This differs from 
(fT4|) at the level of the class of perturbations to be con- 
sidered. We shall come back to the connection between 
the MRS and EHT theories in Sec. UTTCl 



B. Generalized entropies 

We shall now introduce a reduced variational problem 
equivalent to ([36]) but expressed in terms of a generalized 
entropy S[uj] associated with the coarse-grained flow in- 
stead of a functional S x [p] associated to the full vorticity 
distribution. Initially, we want to determine the vortic- 
ity distribution p*(r, a) that maximizes S x [p] with the 
robust constraints E[uJ] = E, T[uJ] = T and the normal- 
ization condition J pda = 1. To solve this maximization 
problem ([36]) . we can proceed in two steps [56|. 

(i) First step: We first determine the distribution 
pi(r, a) that maximizes S x [p] with the constraints E, T, 
J pda = 1 and a fixed vorticity profile uJ(r) = J pa da. 
Since the specification of cJ(r) determines T and E, this 
is equivalent to maximizing S x [p] with the constraints 
J pda = 1 and J pa da = cJ(r). Writing the first order 
variations as 



&(r)5 / pada dr 



C(r)S y pda 



dr = 0, 



where Z(r) and <3>(r) are determined by the contraints 
J pda = 1 and uj = J pa da leading to 



Z(r) = J X (v)e- a * {r) da = x(<S>), 



(41) 



^ J X (cr)ae-^da = -(In *)'(*). (42) 



Z(r) 



J(r) 



Equation ([42]) relates to the vorticity profile uj and Eq. 
([4T]) determines Z. The critical point ([40]) is a maxi- 



mum of S x with the above-mentioned constraints since 
5 2 S X = — J ^^-drda < (the constraints are linear in p 
so their second variations vanish). This gives a distribu- 
tion pi\uj(r),a] depending on u;(r) and a. Substituting 
this distribution in the functional S x [p] , we obtain a func- 
tional S[uJ] = S x [pi] of the vorticity uJ alone. Using Eqs. 
(135]) and (|40|). it is given by 



S[uJ] = / pi(cr$ + lnx)drdcr = / (uJ$ + In dr. 



Therefore, can be written 

5[5J] = - y C(uJ) dr, 



with 



C(uj) 



;<£ - lnx($). 



(43) 



(44) 



(45) 



Now, $(r) is related to cj(r) by Eq. ([42]). This implies 
that 



so that 



C(ZJ) 



-[(lnx)']- 1 (-^), 



(46) 



(47) 



This can be written equivalently 



C(5J) 



/a; 
F~ 1 (x)dx (48) 



where 3>(r) and £(r) are Lagrange multipliers, we obtain 



where the function F is defined by Eq. ([TO]) . Note that 
the function C is convex, i.e. C" > 0. Equation ([44]) with 
([47]) is the entropy of the coarse-grained vorticity. It is 
completely specified by the prior x{°~)- Since the func- 
tion C can take several forms depending on the prior, S 
is sometimes called a generalized entropy [14l4l6|. This 
model of 2D turbulence is therefore an interesting physi- 
(39)cal example where generalized forms of entropy can arise 

Before going further, let us establish some useful iden- 
tities. From Eqs. ([40]) - ([42]) . we easily obtain 



Pi(r,cr) 



Z(r) 



x(o-)e' 



-cr<£(r) 



(40) 



u) 2 ($) = -u/(<&). 



(49) 
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On the other hand, taking the derivative of Eq. ([46]) , we 
have 



1 



dUJ u'(<f>)' 
Combining Eqs. (@9j) and ([50]) . we get [l5|: 

1 



c"{wy 



(50) 



(51) 



Second step: we now have to determine the vortic- 
ity field cJ*(r) that maximizes Sp] with the constraints 
Ep] = E and Yp] = Y. We thus consider the maximiza- 
tion problem 

max {Sp] | Yp] = T, Ep) = E}. (52) 

The critical points of Sp] at fixed E and Y satisfy the 
variational principle 

SS - (3SE - aSY = 0, (53) 

where /3 and a are Lagrange multipliers. This yields 

C'p) = -fii/j-a. (54) 

Using Eq. ([48]) this is equivalent to uJ = + a) and 

we recover the coarse-grained vorticity ([9]) deduced from 
the Gibbs state (|8j). Differentiating the previous relation, 
we note that 



aj'(V>) 



13 



(55) 



According to Eq. ([46]) . we also note that Eq. ([54]) is 
equivalent to 

$ = + a, (56) 
at equilibrium. Then, the identity ([49]) becomes 



^2 --u'{i>)i 



(57) 



returning Eq. ([T2]) . Comparing Eqs. ([55]) and ([57]) . we 
recover Eq. ([5T]) . On the other hand, a critical point of 
32]) is a maximum of S at fixed E and Y iff 1231 : 



£ 2 J = ~\J C"p)(SuJ) 2 dr-^(3 J SuS^dr < 0, (58) 

for all variations SuJ that conserve circulation and energy 
at first order. 

(Hi) Conclusion: Finally, the solution p*(r, a) of ([36]) 
is given by Eq. ([^0]) where cj*(r) is the solution of ([52]) . 
Therefore, p*(r, a) = pi [cJ*(r), a] maximizes 5 x [p] at 
fixed T and normalization iff cJ*(r) maximizes Sp] 
at fixed E and T. Therefore, ([36]) and ([52]) are equiva- 
lent but ([52]) is simpler to study because it is expressed 
in terms of the vorticity field cJ(r) instead of the full vor- 
ticity distribution p(r, a). The equivalence between the 



stability conditions ([38]) and ([58]) is shown explicitly in 
Appendix [B] 

In conclusion, in the EHT approach, the statistical 
equilibrium state p(r, <j) maximizes a relative entropy 
S x [p] at fixed circulation T, energy E and normalization 
condition if and only if the equilibrium coarse-grained 
field cJ(r) maximizes a generalized entropy Sp] (deter- 
mined by the prior) at fixed circulation T and energy E. 
We have the equivalence 



(1361) & ([52]). 



(59) 



This provides a condition of thermodynamical stability 
in the EHT sense. On the other hand, Ellis-Haven- 
Turkington [12| have shown that the maximization prob- 
lem ([52]) also provides a refined condition of nonlinear 
dynamical stability with respect to the 2D Euler equa- 
tion (see [23] for further discussion). Therefore, a EHT 
statistical equilibrium state is both thermodynamically 
stable (with respect to variations of the fine-grained vor- 
ticity distribution Sp) and nonlinearly dynamically stable 
(with respect to variations of the coarse-grained vorticity 
field SuJ). 



C. Another interpretation of the EHT approach 

As noted by Bouchet (22[, and further discussed by 
Chavanis [23], there exists another interpretation of the 
EHT approach. As we have already indicated, the EHT 
approach can be interpreted as a grand microcanonical 
version of the MRS theory in which the Lagrange mul- 
tipliers a n> i are fixed instead of the constraints 
[15j . Therefore, the MRS theory is associated to the mi- 
crocanonical ensemble while the EHT approach is asso- 
ciated to the grand microcanonical ensemble [23|. Now 
it is well-known in statistical mechanics that a solution 
of a maximization problem is always solution of a more 
constrained dual maximization problem j44j|. In particu- 
lar, grand microcanonical stability implies microcanoni- 
cal stability. Therefore, the EHT condition of thermody- 
namical stability provides a sufficient (but not necessary) 
condition of MRS thermodynamical stability. We have 
the implication 



(60) 



This implication can be directly obtained from the stabil- 
ity conditions ([T4]) and ([38]) . Indeed, if inequality ([38]) is 
satisfied for all perturbations that conserve circulation, 
energy and normalization, then it is satisfied a fortiori 
for perturbations that conserve circulation, energy, nor- 
malization and all the Casimirs, so that ([14]) is fulfilled. 
Therefore, an EHT equilibrium is always a MRS equilib- 
rium but the converse is wrong because some constraints 
have been treated canonically. This is related to the 
notion of ensemble inequivalence in thermodynamics for 
systems with long-range interactions (43-[46j. There can 
exist states that solve the maximization problem (|6]) al- 
though they do not solve ([36]) . Such states cannot be 
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reached by a grand microcanonical description. In that 
case, we have ensemble inequivalence. Therefore, an in- 
terpretation of the optimization problem ([36]) is that it 
provides a sufficient condition of MRS thermodynamical 
stability. 

On the other hand, the coarse-grained vorticity field 
associated to a maximum of S[p] at fixed circulation, 
energy, Casimirs and normalization is always a critical 
point of S[uj] at fixed circulation and energy. However, it 
is not necessarily a maximum of S[uJ] at fixed circulation 
and energy since ([52]) is not equivalent to (|6j). Accord- 
ing to ([59]) , the optimization problems ([52]) and ([36]) are 
equivalent so we have the implications 



man & © ©. 



(61) 



Therefore, a maximum of S[uJ] at fixed energy and cir- 
culation is a MRS equilibrium state, but the reciprocal 
is wrong in case of ensemble inequivalence. For example, 
if the MRS equilibrium vorticity distribution is Gaussian 
(which corresponds to specific initial conditions), the gen- 
eralized entropy S[uJ] is proportional to minus the coarse- 
grained enstrophy T^' 9 ' = JuJ 2 dr (see Sec. [V)). There- 
fore, a minimum of coarse-grained enstrophy at fixed en- 
ergy and circulation is a MRS equilibrium state, but the 
reciprocal is wrong in case of ensemble inequivalence. 

Remark 1: In case of equivalence between microcanon- 
ical and grand microcanonical ensembles, these results 
justify a "generalized selective decay principle" [23|. In- 
deed, in that case, the equilibrium coarse-grained vortic- 
ity field maximizes a generalized entropy S[oJ] (or mini- 
mizes the functional — S) at fixed circulation and energy. 
For a Gaussian equilibrium vorticity distribution, this 
justifies a minimum coarse-grained enstrophy principle 
through statistical mechanics. In the present case, the 
increase of generalized ent ropy is due to coarse-graining: 
the microscopic Casimirs S[lj] calculated with the fine- 
grained vorticity are conserved while the macroscopic 
Casimirs S\uJ] calculated with the coarse-grained vortic- 
ity increase. By constrast, the energy E[uJ] and the cir- 
culation T[uJ] calculated with the coarse-grained vorticity 
remain approximately conserved. 

Remark 2: Since (|6)) is not equivalent to ([52]) , a MRS 
statistical equilibrium state does not necessarily satisfy 



the condition of refined dynamical stability ([52]) given by 
Ellis-Haven- Turkington [12j. However, it can be shown 
that a MRS statistical equilibrium state is always non- 
linearly dynamically stable with respect to the 2D Euler 
equations as a consequence of the Kelvin-Arnol'd theo- 
rem which provides an even more refined condition of 
nonlinear stability than the EHT criterion (see [23| for 
details). 



IV. RELAXATION EQUATIONS WITH A 
PRIOR VORTICITY DISTRIBUTION 

Let us now derive some relaxation equations associated 
with the EHT approach. These relaxation equations will 



be compared to those associated with the MRS theory. 



A. A first type of relaxation equations 

Chavanis [lil-HiSj has proposed a relaxation equation 
solving the maximization problem ([36]) . This can serve 
as a numerical algorithm to compute maximum entropy 
states with appropriate constraints. The idea is to use 
the two-steps method presented in Sec. IIIIBI The dis- 
cussion is here slightly improved. We assume that, at any 
time of the evolution, the vorticity distribution p(r, a, t) 
maximizes the relative entropy ([35]) with the constraint 
on mean vorticity cJ(r,t) = f pa da and normalization 
J pda = 1. This leads to the time dependent distribu- 
tion 



-<r$(r,t) 



Z{r,t)' 

where Z(r,t) and <3>(r, t) are determined by 



(62) 



Z(r,t) 



x{a)e -«Hr,t) da = ^ 



(63) 



"(r,t) = / xi^ae-^Ua = -(ln*)'(*). 

(64) 

Now, according to (jUJ), we know that the equilibrium 
vorticity field tU(r) maximizes the generalized entropy 



with 



S[u 



C(ZJ) = 



C(u) dr, 



[(InxYrH-x^x, 



(65) 



(66) 



at fixed circulation and energy. We can obtain a re- 
laxation equation for cJ(r, t) solving this maximization 
problem by using a generalized Maximum Entropy Pro- 
duction Principle [14Ml6j. We assume that the coarse- 
grained vorticity evolves in time so as to maximize the 
rate of (generalized) entropy production S[uJ] (fixed by 
the prior) at fixed circulation and energy. This leads to 
a generalized Fokker-Planck equation of the form 



duo 
~dt 



u • VcJ = V • < D 



with 



j DVuJ ■ V?/><fr 



(67) 



(68) 



where (3{t) is a Lagrange multiplier enforcing the energy 
constraint E = at any time and D(r, t) is the diffusion 
coefficient. The boundary conditions are 



Vw ■ 



C"(U) 



n 



0, 



(69) 
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where n is a unit vector normal to the boundary, in or- 
der to guarantee the conservation of circulation. Easy 
calculations give the generalized i^-theorem 



J DC"(W) (vuJ 



dr>0, (70) 



C"(5J) 



provided that D(r,t) > 0. By construction, the relaxed 
vorticity field uJ(r) solves the maximization problem ([52]) . 
Then, the corresponding distribution ([62]) - ([M]) solves the 
maximization problem ([36]) . Therefore, these relaxation 
equations tend to the statistical equilibrium state corre- 
sponding to the EHT approach. The diffusion coefficient 
is not given by the MEPP, but it can be estimated by 
Eq. ([18]) . Using Eqs. ([62])- ([64]) and repeating the steps 
([45]) - ([5T]) . we establish that at any time(57l]: 



0J2 



c"(ujy 



(71) 



Plugging this result in Eq. ([18]) . we obtain the expression 
of the diffusion coefficient 



D(r,t) 



(72) 



Finally, using arguments similar to those of Sec. Ill Dl 
the equation for the velocity field is 



du 

~dt 



(U • V)U : 



: Vp + D Au 



C"(57) 



At equilibrium, we get 



Au 



C"(uJ) 



— u, 



• (73) 



(74) 



which can be directly derived from the Gibbs state ([8]) 
using Au = z x VcJ, VcJ uj' (t/j)\7tp and Eq. ([55]) . 

Therefore, the system of equations ([62]) - ([72]) is com- 
pletely closed when the prior vorticity distribution x( a ) 
is given [5^. Note that they determine not only the evo- 
lution of the mean flow cJ(r, t) through Eqs. ([67]) . ([68]) 
and ([72]) but also the evolution of the full vorticity dis- 
tribution p(r,a,t) through Eqs. ([62]), ([63]) and dM}. The 
case of a Gaussian prior will be discussed specifically in 
Sec. [V] However, we stress that the relaxation equations 
([62]) - ([72]) . which are relatively easy to solve numerically, 
can be used to study situations going beyond the Gaus- 
sian approximation. Some examples showing the con- 
struction of the generalized entropy S\uJ] from the prior 
x(&) and giving the corresponding equilibrium states and 
the corresponding relaxation equations are discussed in 
[lil-HiSj . Apart from their potential interest in 2D turbu- 
lence, this leads to interesting classes of nonlinear mean 
field Fokker-Planck equations [18| . 

It should be emphasized that the relaxation equations 
associated with the maximization problem ([52]) are not 



unique. For example, another type of relaxation equa- 
tions (see Appendix [D| solving the maximization prob- 
lem ([52]) is given by [23[ : 



duo „ 
— - + u • VcJ 

at 



-D(C(u) + P(t)4> + a(t)), (75) 



where D is a positive coefficient and the Lagrange mul- 
tipliers j3(t) and a(t) evolve according to 



(DC(uj)^) + p(t){D^ 2 ) + a(t)W> = 0, 



(76) 



(DC (W)) + P{t)m) + a(t)(D) = 0, (77) 

so as to conserve the energy and the circulation (the 
brackets represent the domain average (X) = J Xdr). 
The boundary conditions are 



C'(ZJ) = -a(t), 



(78) 



on the boundary, in order to be consistent with the equi- 
librium state where the r.h.s. of Eq. ([75]) is equal to zero 
on the whole domain (recall that ip = on the boundary) . 
Easy calculations lead to the generalized i^-theorem 



S= D(C(u) + f3(t)ilj + a(t)) 2 dr>0. 



(79) 



The relaxed vorticity field cJ(r) solves the maximization 
problem ([52]) . Then, the corresponding distribution ([62]) - 
([64]) solves the maximization problem ([36]) . 

Remark: Since the maximization problem ([52]) pro- 
vides a refined criterion of nonlinear dynamical stability 
for the 2D Euler equations [12[ , the relaxation equations 
of this section can also be used as numerical algorithms to 
construct nonlinearly dynamically stable steady states of 
the 2D Euler equation independently from the statistical 
mechanics interpretation [231 ] . 



B. A new type of relaxation equations 

In the previous approach, the vorticity distribution is 
assumed to have the form ([62]) - ([64]) at each time. This 
assumption is consistent with the two-steps method de- 
veloped in Sec. IIII Bl to show the equivalence between 
the maximization problems ([36]) and ([52]) . We shall now 
introduce a new type of relaxation equations in which the 
form of the vorticity distribution changes with time [23| . 
These relaxation equations are associated with the basic 
maximization problem ([36]) . 

We write the equation for the evolution of the vorticity 
distribution p(r, cr, t) in the form 



dp „ dJ 



(80) 



where J(r, a, t) is a current acting in the space of vor- 
ticity levels a. It will be determined by a systematic 
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procedure. We note that this form assures the conserva- 
tion of the local normalization J pda = 1 provided that 
J — >> for a — » ±oo. We also emphasize that the to- 
tal areas of the vorticity levels 7(0-, t) = J pdr are not 
conserved by the relaxation equations (|80j) . This is be- 
cause, in the EHT approach, the Casimirs are not con- 
served. This differs from the relaxation equations ([T5]) 
associated with the MRS approach where the left hand 
side is of the form —V • J where J(r, cr, t) is a current 
acting in position space. In the MRS approach, the 
Casimirs, or equivalently the total areas of the vortic- 
ity levels 7(0, t) = j pdr, are conserved by the relaxation 
equations. This is not the case in the EHT approach. 
Therefore, we see from the start that the structure of 
the relaxation equations will be very different in the two 
approaches. 

Multiplying Eq. ([80]) by a and integrating on the vor- 
ticity levels, we obtain an equation for the coarse-grained 
vorticity 



duo 
~dt 



u • Vco> 



/ 



Jda, 



(81) 



where we have used an integration by parts to get the 
r.h.s. We can also derive an equation giving the evolution 
of the local centered variance 002- Multiplying Eq. ([80]) 
by a 2 and integrating on a, we obtain 



duo 2 



u • Vuj 2 = 2 Jcr da, 



(82) 



where we have used an integration by parts to get the 
r.h.s. On the other hand, multiplying Eq. ([8T]) by ZU, we 

get 



u • Vuj 2 = 2uu I J da. 



duo 2 
dt 

Subtracting these two equations, we find that 



(83) 



^ ± u • Voj 2 = 2 J J {a - uJ)da. (84) 

Let us now consider the constraints that the relaxation 
equations must satisfy. The conservation of the circula- 
tion implies 



r = J Jdrda = 0, (85) 

and the conservation of the energy implies 

E = J Jtpdrda = 0. (86) 

On the other hand, the rate of production of relative 
entropy S x is given by 



5 V 





p 







drda. 



(87) 



We shall determine the optimal current J by maximizing 
the rate of relative entropy production at fixed circulation 
and energy. We also introduce the physical constraint 
J J 2 /(2p)do~ < C that puts a bound on the current. We 
write the variational problem as 

5S X - P(t)5E - a(t)5T - J -^5 (i^ dc ^) dr = °> 

where a(t), j3(t) and D(y, t) are Lagrange multipliers that 
assure the conservation of circulation and energy at any 
time. Performing the variations, we find that the optimal 
current is 



J 



-Dp 



\da 



In 



Lx(<x)J 



+ 0{t)1> + a{t)\, (89) 



or equivalently 
dp 



J 



da 



p(\nxy(<T)+l3(t)p4> + a(t)p\- (90) 



Therefore, the relaxation equation for the vorticity dis- 
tribution p(r, a, t) takes the form 



dp 



d_ 

da 



P (ln x)'{p) + P(t)p4> + <*{t)p 



.(91) 



Using Eq. (j8T]h we obtain the relaxation equation for the 
coarse-grained vorticity 



— - ± u • VcJ = D 

at 



p(]nx)'((r)d(T - P(t)i) - a(t) 



(92) 



In general this equation is not closed as it depends on 
the full vorticity distribution p(r,cr, t). The Lagrange 
multipliers are determined by Eqs. ([85]) . ([86]) and ([90]) 
yielding 

j (Dp)OnxY((T)dcT - mm) - a(t)(D) = 0, (93) 
J (D p^) (In X y(o-)da - f3(t)m 2 ) - a(t)(D^) = 0. (94) 



Substituting Eq. ([89]) in Eq. ([87]) and using the con- 
straints ([85]) - ([86]) . it is easy to establish the H-theorem 



J Dp 



drda > 0, 



(95) 



provided that D > 0. On the other hand, a stationary 
solution corresponds to J = yielding 



do- [ x(o-) 



± ± a = 0. 



(96) 
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After integration, we recover the Gibbs state (j8j). Equa- 
tion ([9T]) with the constraints ([93]) - ([M]) has the following 
properties: (i) T, E and J pda = 1 are conserved, (ii) 
S x > 0. (hi) S x = p(r, a) is the Gibbs state (JSJ) 
<^> 9tp = 0. (iv) p(r, a) is a steady state of Eqs. ([9T]) . 
([93]) . ([94]) iff it is a critical point of *S X at fixed T, E 1 and 
normalization, (v) a steady state of Eqs. ([9T]) . ([93]) . ([94]) 
is linearly stable iff it is a (local) maximum of S x at fixed 
T, E and normalization. By Lyapunov's direct method, 
we know that if S x is bounded from above, Eqs. ([91]) . 
([93]) . ([94]) will relax towards a (local) maximum of S x 
at fixed T, £" and normalization (if several local maxima 
exist, the choice of the maximum will depend on a notion 
of basin of attraction). Therefore, a stable steady state 
of Eqs. ([9T]) . ([93]) . ([94]) solves the maximization problem 
([36]) for a given prior and a given circulation and energy 
specified by the initial condition. 

We can easily derive a hierarchy of equations for the 
local moments of the vorticity. Multiplying Eq. ([80]) by 
a n and integrating on a, we get 



duo n 



• u • Vuo n 



Ja n ~ Y da. 



(97) 



Inserting the expression ([90]) of the current, we obtain 



dt 



+ u • Vu n = Dn(n - l)uj n - 2 



If we recall that 



lnx(cr) = -J2 



a k a 



(98) 



(99) 



k>l 



we can rewrite the foregoing equation in the form 



duo n 



u • Vuj u = Dn(n - l)uj n 



Dnl > a k kuj k + n ~ 2 + pit)^- 1 ^ + a^cj 71 ' 1 }. (100) 



From this general expression, we see that the hierarchy of 
equations is closed only if a k = for k > 2 that is to say 
for a Gaussian prior (see next section). Otherwise, the 
equation for the moment of order n requires the knowl- 
edge of the moment of order n + (k — 2) > n and some 
closure approximations must be introduced. This can be 
an interesting mathematical problem but it will not be 
considered in this paper. 

Remark: According to Eq. ([60]) . the relaxation equa- 
tions presented in this section also solve the dual maxi- 
mization problem (j6j) for the corresponding values of the 
Casimirs. Therefore, according to the interpretation of 
Sec. IIII CI they can be used as numerical algorithms to 
construct a subclass of MRS statistical equilibria. 



V. THE CASE OF A GAUSSIAN PRIOR 

In this section, we consider the particular case of a 
Gaussian prior vorticity distribution and make the con- 
nection with minimum enstrophy states. 

A. Equilibrium states 

We assume a Gaussian prior of the form 



i 



(101) 



where ^2 is a constant. The corresponding Gibbs state 
(El) is 



1 1 



(102) 



The vorticity and the local centered variance are given 

by 



and 



5J(r) = -n 2 W(r)+a), 



w 2 (r) = ft 2 - 



(103) 



(104) 



Therefore, in the case of a Gaussian prior, the ZZJ — ip 
relationship is linear and the local centered variance is 
uniform. The Gibbs state (|102p can be rewritten 



1 

pM = vm e 2(12 



(105) 



At equilibrium, substituting Eq. (|105|) in Eqs. ([2]) and 
{35}, we get S[p eq ) = \ ln(lf * - T c 2 9 ') and S x [p eq ] = 
~ 2^1 ^2 9 ( U P t° additive constants). On the other hand, 
the generalized entropy defined by Eqs. (|45j) and ([47]) 
can be conveniently calculated from Eqs. ([54]) and (j!03p 
yielding C'(uJ) = hence 



5: 



2a 



-dr 



2Q 2 2 



(106) 



Therefore, for a Gaussian prior, the generalized entropy is 
proportional to minus the macroscopic enstrophy Y c 2 9 ' = 
f uj 2 dr. In that case, according to Eq. ([59]) . the maxi- 
mization of the relative entropy S x at fixed energy and 
circulation (EHT thermodynamical stability) is equiva- 
lent to the minimization of the macroscopic enstrophy 
T^' 9 ' at fixed energy and circulation, i.e. 

min {r c 2 g -p} 1 rp] = r, Ep] = E}. (107) 



According to the interpretation given in Sec. IIII CI this 
minimization problem also provides a sufficient condition 
of MRS thermodynamical stability. 
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Remark: Writing the variational problem in the form 
([53]) , the critical points of 5, given by Eq. (jl06p , at fixed 
energy and circulation are given by Eq. (j!03|) . Further- 
more, they are maxima of S at fixed energy and circula- 
tion iff 

5 2 J = --^- [(SuJ) 2 dr-lp [ SuSi/jdr < 0, (108) 

for all perturbations JcJ that conserve energy and circu- 
lation at first order. 



B. A first type of relaxation equations 



also that the equation for the velocity field ([73]) takes the 
form 



du 



(u- V)u 



1 



dt * p 
At equilibrium, we get 

Au = f30 2 u. 



Vp + D (Au - /3(t)0 2 u). (116) 



(117) 



On the other hand, using C'(uJ) = ZJ/O2 according to 
Eq. (|106p . the alternative relaxation equations ([75]) - ([77]) 
become (we here assume D constant to slightly simplify 
the expressions) 



We first consider the relaxation equations of Sec. IIV Al 
In this approach, for a Gaussian prior, the vorticity dis- 
tribution is at any time given by 



du „ 
— - + u • VZJ 

dt 



-D 



-^uJ + P(t)tl> + a(t) 



(118) 



1 



(<x-ZJ(r,t))^ 



v/2^ 



(109) 



Therefore, the vorticity distribution is Gaussian and the 
local centered variance cj2(r,t) = ^2 is uniform and con- 
stant in time. Using C"(uJ) = 1/0,2 according to Eq. 
(|106|h the evolution of the vorticity is given by 



9oj „ „ 
— - + u • VcJ = V • 
at 



D 



]}■ 



/ DVuJ ■ Vipdr 
7 I>fi 2 (V^) 2 dr' 



(110) 



(111) 



a(t) 



1 T{ip) - 2AE 



1 r(^ 2 ) - 2E{tp) 



(119) 



(120) 



where A is the domain area. The i^-theorem ([79]) can be 
written 



1 

20 2 



I 



I ^- + (3{t)ip + a(t) 



(121) 



D = Ke 2 n 2 . 



(112) 



C. A new type of relaxation equations 



Since the diffusion coefficient D is constant, Eqs. (|110|) 
and (jllip can also be written 



duo 
~dt 



u • VcJ = D(AcJ - /3(t)fi 2 aJ), (113) 



fuo 2 dr _ -T c 2 g -(t) _ S(t) 



2 J(Vij) 2 dr 20 2 E 



E ' 



(114) 



where we have used an integration by parts to simplify 
the last equation. Interestingly, for this Gaussian model, 
we see that S{t) = (3(t)E out-of-equilibrium. On the 
other hand, the H theorem can be written 

^ = 'wfi 9 ' = I £ (VZZJ + w^ v ^ 2 dv ^ °- 

(115) 

Therefore, the relaxation equations increase the general- 
ized entropy (or decrease the coarse-grained enstrophy) 
at fixed energy and circulation until the maximum en- 
tropy (or minimum enstrophy) state is reached. Note 



Let us now consider the new relaxation equations of 
Sec. IIVBI Since (lnx)'(c") = — &/O2 for the Gaussian 
prior (jlOlj) . the current ([90]) becomes 



J 



-D 



1 



dp 
d^^O: 



pa + f3(t)pi/j + a(t)p 



(122) 



and the evolution of the vorticity distribution is given by 
a relaxation equation of the form 



d_ 

~ da 



u- Vp 



(123) 



We note that the form of the distribution p(r, cr, t) 
changes with time. The vorticity distribution is Gaus- 
sian only at equilibrium. The relaxation equation for the 
coarse-grained vorticity ([92]) becomes 



duo ^ 
— - + u • VcU 

dt 



-D 



-U; + p(t)i/> + a{t) 

\L2 



(124) 
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and the equations (|93|) - (|94|) for the Lagrange multipliers 
reduce to (we here assume D constant to slightly simplify 
the expressions) 



/3(f) 



a(t) 



1 T(ip) - 2AE 
ThA{^) - <V>> 2 



i r(y, 2 ) - 2E{*i>) 

Q 2 A(ip 2 } - (4>) 2 ' 



(125) 



(126) 



We emphasize that, for a Gaussian prior, the equation for 
the coarse-grained vorticity ([124)) is closed. Furthermore, 
it coincides with Eq. ([118)) . Using Eq. ([Hi]) and the 
expression (j!22p of the current, we find that the equation 
for the local centered variance is 



duj 2 



u • Vo; 2 = 2D[1 



0J2 



(127) 



This is also a closed equation. Using Eq. (jlOQj) . the relax- 
ation equations for all the local moments of the vorticity 
are 



■ u • Vuj n = Dn(n - l)u n 



-Dn 



\L2 



(128) 



Finally, for a Gaussian prior vorticity distribution, we 
can explicitly check that the generalized entropy ([106)) 
increases monotonically for a dynamics of the form ([123j) . 
Indeed, since Eqs. ([124)) and ([118)) coincide, we directly 
obtain the i^-theorem ([121)) . We have not been able to 
prove that the generalized entropies S[uJ] increase mono- 
tonically for a dynamics of the form ([9T]) in the case of a 
prior that in not Gaussian. This could be an interesting 
mathematical problem to consider. 

Remark: Note that D plays the role of a diffusion co- 
efficient in the space of vorticity levels [see Eq. ([123)) ]. 
However, in the vorticity equation ([124)) , the quantity 
D/Q.2 ~ l/trelax plays the role of an inverse relaxation 
time. 



VI. NUMERICAL SIMULATIONS 

In this section, we generalize the relaxation equations 
derived in Sec. IV CI so as to account for the presence 
of a bottom topography h(x, y) and we solve them nu- 
merically to illustrate the nature of phase transitions in 
geophysical flows. 



A. The quasigeostrophic equations 




FIG. 1: Continental slope: (a) topography h(x,y); (b) to 
(d) Equilibrium PV field for (b) /5 = 1000, (c) f3 = 10, (d) 
(3 = -10. 



R — >• +oo. It is described by the quasigeostrophic (QG) 
equations 

^ + u-V<? = 0, q = -A^ + /z, u = -z x V^, (129) 

where q is the potential vorticity and uz = V x u is 
the vorticity satisfying uj = —Aifj. The QG equations 
conserve the energy 

E=\J(q- h)i>dv, (130) 



and the Casimirs 



If = J f(q)dr, 



(131) 




Let us consider a 2D incompressible flow over a to- 
pography h(x,y) in the limit of infinite Rossby radius 



FIG. 2: Seamounts: (a) topography h(x,y); (b) to (d) Equi- 
librium PV field for (b) /5 = 1000, (c) = 10, (d) f3 = -10. 
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FIG. 3: Ridge: (a) topography h(x,y)\ (b) to (d) Equilibrium 
PV field for (b) {3 = 1000, (c) f3 = 10, (d) f3 = -10. 

where / is an arbitrary function. We consider a general- 
ized entropy of the form 

s '-wJ fit - (132) 

corresponding to a Gaussian prior (see Sec. |Vj). As shown 
in [2l| , the critical points of entropy (j!32p at fixed energy 
(|130p and potential circulation T = j qdr are solutions of 
the differential equation 

- + W = T + f3(<iP) - ft, (133) 

with tp = on the domain boundary. 

We thereafter illustrate the fact that the relaxation 
equations derived in Sec. IV CI can be used as numerical 
algorithms to determine maximum entropy states. We 
consider three topographies h(x , y) similar to those in- 
troduced by Wang & Vallis [47j: continental slope (see 
Fig. Ufa)), seamounts (Fig. E{a)) and ridge (Fig. EJa)). 

B. Equilibrium states 

As a first step, we show in Figs. [TJ El and El the coarse- 
grained potential vorticity at equilibrium for T = and 
different values of (3 in a square domain. These vortic- 
ity fields have been obtained by solving Eq. (|133j) with 
the method given in [21[. As pointed out later, all these 
states are stable (they correspond to maximum entropy 
states for the corresponding values of the energy). Two 
remarkable trends can be observed: (i) for very large val- 
ues of /3, the potential vorticity has the tendency to align 
with the topography. This corresponds to generalized 
Fofonoff flows (<ip ~ -ft//3, q ~ ft, u ~ ±z x Vft) HH; (ii) 
the PV fields for f3 > and f3 < are of opposite signs. 

More quantitatively, we plot in Fig. 2] the curve 
/3(l/i£)[59j for the "seamounts" topography in a square 
domain with T = 0. For small energies, Eq. (|133p admits 



only one solution, whereas for larger values of E it admits 
an infinite number of solutions, i.e. there exists an infi- 
nite number of critical points of entropy at fixed energy 
and circulation. In the latter case, the maximum entropy 
state is the one with the highest f3 pjHIlj . For small en- 
ergies, we have Fofonoff flows (j3 > 0), for intermediate 
energies we have reversed Fofonoff flows (j3 < 0), for large 
energies we have a direct monopole and for E — >> +00 we 
have a monopole rotating in either direction with inverse 
temperature A detailed description of this caloric 

curve can be found in [20|,l21[. Since r ^ T* (where T* is 

defined in [2lj|). there is no plateau at f3^\ Therefore, the 
monopole is obtained smoothly from the Fofonoff flows 
in the limit E —> +oo. In particular, in the present sit- 
uation, the caloric curve /3(E) does not display a second 
order phase transition. 

On the other hand, Chavanis & Sommeria [l9| have 
demonstrated that 2D Euler flows characterized by a lin- 
ear uj — i/j relationship experience geometry induced phase 
transitions between a monopole and a dipole when the 
domain becomes sufficiently elongated. In the case of a 
rectangular domain, this transition occurs at a critical 
aspect ratio r c = 1.12 (the maximum entropy state is the 
monopole for l/r c < r < r c and the dipole for r > r c 
or r < 1/t c ). As shown by Venaille & Bouchet [20| and 
Naso et al [2l|, this property persists in the case of QG 
flows with a topography. To illustrate this result, we plot 
in Figs. E] and [6] the curve f3(l/E) for the "seamounts" 
topography in a rectangular domain with T = 0. We 
have considered rectangular domains with aspect ratios 
r = 2 > r c (horizontal) and r = 1/2 < l/r c (vertical) 
in order to emphasize the difference with respect to the 
case of a square domain (l/r c < r = 1 < r c ). For small 
energies, we have Fofonoff flows (j3 > 0), for intermedi- 
ate energies we have reversed Fofonoff flows ((3 < 0), for 
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1/E 

FIG. 4: Relationship between /3 and 1/E in a square domain 
with the topography of Fig. EI a) and r = 0. The q density is 
plotted for several values of f3 (increasing values from blue to 
red). 
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C. Relaxation towards the maximum entropy state 



FIG. 5: Relationship between f3 and 1/E in a rectangular 
domain with the topography of Fig. [2ja) and r = (r = 2). 




2000 4000 6000 8000 10000 12000 



FIG. 6: Relationship between f3 and 1/E in a rectangular 
domain with the topography of Fig.[2ja) and r = (r = 1/2). 



large energies we have a direct dipole and for E — » +00 
we have a dipole rotating in either direction with inverse 
temperature ^21 (horizontal) or j3\2 (vertical). A detailed 
description of these caloric curves can be found in (201121). 
Since the eigenmodes (2,1) and (1,2) are not orthogo- 
nal to the topography, there is no plateau at fci or f3\2. 
Therefore, the dipole is obtained smoothly from the Fo- 
fonoff flows in the limit E — »• +00. In particular, in the 
present situation, the caloric curve /3(E) does not display 
a second order phase transition. 

As a final remark, we note that low energy states are 
strongly influenced by the topography (and only weakly 
by the domain geometry) as in the FofonofT [48j study 
whereas high energy states are strongly influenced by the 
domain geometry (and only weakly by the topography) 
as in the Chavanis-Sommeria [l9| study. 



Following the same approach as in Sec. II V Bl and con- 
sidering a Gaussian prior like in Sec. IV CI one can derive 
the following relaxation equations for the coarse-grained 
potential vorticity and its local centered variance 



dq 
dt 



u-Vq 



dt 



-D 



_Q_ 

Q2 



U • V<?2 



- /3(t)<i/) + a(t) 
2D ( 1 



q2_ 
Q2 



The conservation of E and T implies 
1 T(^) - A(2E-{ 



f3(t) = 



a(t) 



Q2 A{^) - ' 

1 r<^ 2 >- MpE + tyh)) 



Q2 



(134) 
(135) 

(136) 
(137) 



We thereafter integrate numerically Eqs. (|134p ~ (|137p 
with the following boundary conditions 



^\dD = 0, 

q\dD -Q2<*(t), 
q2\dD Q2, 



(138) 
(139) 
(140) 



where dD is the domain boundary. In the simulations, 
all the quantities are normalized by the lengthscale A 1 / 2 

— 1/2 

and by the timescale Q 2 (this amounts to taking A = 
Q2 = 1 in the foregoing equations). Furthermore, the 
potential circulation T is taken equal to and the energy 
is normalized by 2/b 2 . 

We first illustrate the convergence of these relax- 
ation equations towards the maximum entropy state, 
integrating them with initial condition: q(x,y,t = 
0) = qo sin(47rx) sin(47n/) and q2(x,y,t = 0) = 
sin 2 (27rx) sin 2 (27n/) + Q2- Varying the value of the pa- 
rameter qo enables us to generate initial conditions of 
different energy E. The dynamics of q and q2 are not 
coupled (except through the advective term). We there- 
fore do not show here the time evolution of q2 that always 
converges to the uniform state q2(x,y) = Q2. 

We plot in Fig. [3 the time evolution of f3(t) and g(r, t) 
for three values of E in a square domain. For each value 
of energy E, the relaxation equations converge towards 
the corresponding maximum entropy state (see plateaus 
in Fig. [7]). The value of f3 in the final state, thereafter 
denoted /3f, is in every case in agreement with the value 
reported in Fig. 3J For large values of 1/E (Fig. Ufa) 
and (b)), the relaxation equations converge towards a Fo- 
fonofT flow which is their unique steady state. A case of 
interest is the one arising for 1/E = 0. In that case, there 
exists an infinite number of steady states and the initial 
condition that we have imposed is an unstable steady 
state. After some time, the system destabilizes and fi- 
nally converges towards the (stable) maximum entropy 
state. In a square domain, the maximum entropy state of 
infinite energy is the monopole for which (3 = (3^. The 
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FIG. 7: Relaxation of the system at fixed E, with the to- 
pography of Fig. [2{a), D = 1 and r = (square domain): 

(a) 1/E « 6620, f3 f « -2.1 (see vertical line in Fig. [4j), 

(b) l/£ « 1196, p f « -34 (see vertical line in Fig. g)), (c) 

1.3.10- 5 , 



limit of infinite energy cannot be reached numerically, 
but we can approach it (see Fig. 0(c)). As expected, for 
very large but finite values of E, the stable steady state 
is the direct monopole which is the natural evolution of 
a Fofonoff flow as energy increases. We plot in Fig. [8] the 
time evolution of j3(t) and g(r, t) in a rectangular domain 
with aspect ratio r = 2 > r c . For a sufficiently large en- 
ergy, the system evolves towards a dipole instead of a 



monopole. 




t 



FIG. 8: Relaxation of the system at fixed E, with the topog- 
raphy of Fig. H D = 1, T = 0, r = 2 and 1/E & 1.6.1CT 5 , 
Pf «#n (see Fig. O. 

We now investigate the influence of the relaxation term 
(r.h.s. in Eq. ([134)) ) on the dynamics. To that purpose, 
we integrate the relaxation equations with different val- 
ues of D, starting from the same initial condition: a 
random field written as the sum of sine functions with 
random phases and amplitudes, and wave numbers rang- 
ing from 1 to 5. As shown in Fig. [91(a), corresponding 
to an energy 1/E = 7200, the system always relaxes to- 
wards the maximum entropy state, with Pf « 1.83. As 
expected, this state is reached at longer times when D 
is decreased. More interestingly, D measures the rela- 
tive importance of the relaxation term (that forces the 
convergence towards the maximum entropy state) and 
of the advection (related to the mixing). Therefore, the 
intermediate vorticity fields are not identical for differ- 
ent values of D. For a smaller D, the advection can 
play a larger role (i.e. the mixing is more efficient) be- 
fore the final state is reached. This behavior is clearly 
seen by comparing Figs. [5fb) and (c). If we use the 
relaxation equations as a parametrization of 2D turbu- 
lence (see, however, the last paragraph of the conclu- 
sion), we see that the influence of the relaxation term is 
to smooth out the small scales without strong influence 
on the large-scale dynamics. This is precisely the role of 
a parametrization since we are in general interested by 
the largest scales, not by the fine structure of the flow. 
Our proposed parametrization rigorously conserves the 
circulation and the energy (contrary to a parametriza- 
tion involving only a turbulent viscosity or an hypervis- 
cosity) and "pushes" the system towards the statistical 
equilibrium state with Gaussian fluctuations. Therefore, 
it allows to simulate 2D flows with less resolution than 
usually required (the simulations have been performed 
on a 128 2 grid) since the small scales have been modelled 
in an optimal manner. It would be interesting, however, 
to compare the efficiency of different par ametrizat ions. 
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FIG. 9: Relaxation of the system at fixed E, with the to- 
pography of Fig. [2{a) and r = (square domain): (a) 
D = 0.1, 0.2, 0.4, 1; (b) D = 0.1; (c) D = 1. For small values 
of D, the mixing is efficient and leads to filament-like struc- 
tures in the transient states. 



VII. CONCLUSION 

In this paper, we have proposed a new class of 
relaxation equations associated with the Ellis-Haven- 
Turkington statistical theory of 2D turbulence where a 
prior vorticity distribution is prescribed instead of the 
Casimir constraints in the Miller- Robert- Sommeria the- 
ory. We have considered specifically the case of a Gaus- 



sian prior associated with minimum enstrophy states and 
we have given a numerical illustration of these relaxation 
equations in connection to Fofonoff flows in oceanic circu- 
lation. We have discussed the connections with previous 
relaxation equations introduced by Robert & Sommeria 
[l3[ and Chavanis [3-16]. These relaxation equations 
can provide efficient numerical algorithms to solve the 
various constrained maximization problems appearing in 
the statistical mechanics of 2D turbulence [23[. This is 
clearly an interest of these equations because it is usu- 
ally difficult to directly solve the Euler-Lagrange equa- 
tions ([T3]) for the critical points of entropy and be sure 
that they are entropy maxima. The relaxation equations 
assure, by construction, that the relaxed state is a maxi- 
mum entropy state with appropriate constraints. In this 
sense, the relaxation equations can provide an alternative 
and complementary method to the numerical algorithm 
of Turkington & Whitaker [49] that has been introduced 
to solve constrained optimization problems. An advan- 
tage of using relaxation equations is that, by incorpo- 
rating a space and time dependent diffusivity related to 
the strength of the fluctuations [3l|, [33| , one can heuris- 
tically account for incomplete relaxation [H, [39| , which 
is not possible with the Turkington- Whitaker algorithm. 
Furthermore, even if these relaxation equations cannot 
be considered as a parametrization of 2D turbulence (see 
below) , they may nevertheless provide an idea of the true 
dynamics towards equilibrium. In that respect, it would 
be interesting to compare them with large eddy simu- 
lations (LES). These relaxation equations also provide 
new classes of partial differential equations which can be 
of interest to mathematicians. In fact, several mathe- 
matical works have started to study these classes of re- 
laxation equations [34]- 36] . The connection to nonlinear 
mean field Fokker-Planck equations and to the Keller- 
Segel model of chemotaxis in biology is also interesting 
to mention in that respect [l8j . If we want to address 
more fundamental issues from first principles, we can try 
to develop kinetic theories such as the quasilinear theory 
of the 2D Euler equation [40j| , the rapid distortion theory 
[50| or the stochastic structural stability theory (5l) . In 
that respect, we would like to briefly comment on the H- 
theorem in 2D turbulence. Ideally, the kinetic equations 
for p(r, a, t) and cJ(r, t) should be derived from first prin- 
ciples, as attempted in [ioj], and the ^-theorem should 
be derived from these kinetic equations. Here, we have 
used the inverse approach: assuming that an i^-theorem 
holds for some functionals, we have constructed relax- 
ation equations that satisfy this i^-theorem while con- 
serving some constraints. This is clearly a purely ad hoc 
procedure that is sufficient to construct numerical algo- 
rithms solving constrained maximization problems, but 
not sufficient to conclude that these relaxation equations 
constitute an accurate parametrization of 2D turbulence. 
Our main goal, here and in [23|, was to provide a relax- 
ation equation for each optimization problem introduced 
in 2D turbulence. Therefore, there exists as many re- 
laxation equations as maximization problems. These re- 
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laxation equations can be used as numerical algorithms 
to solve these optimization problems and therefore de- 
termine dynamically and/or thermodynamically stable 
steady states. This is the main virtue of these relaxation 
equations [23j. 



Appendix A: Maximization of entropy in three steps 

Let us consider the maximization problem (see also 
0): 

max {S[p\ | Tp] = T, Ep] = E, 

p 

rf<V] = r^, J P da = i}, (Ai) 

where S is the MRS entropy ((2]). To solve the maximiza- 
tion problem (|A1[) . we can proceed in three steps. This 
will show the connection with the results of Sec. Ill CI 

First step: we first maximize S[p] at fixed E, T, T^' 9 ', 
normalization and a fixed profile of vorticity cJ(r) = 
J pa da and enstrophy uo 2 (r) = J pa 2 da. Since the spec- 
ification of uo(r) and uo 2 (r) determines E, T and T^', 
this is equivalent to maximizing S[p] at fixed normaliza- 
tion for a fixed profile of vorticity uo(r) = j pa da and 
enstrophy uo 2 (r) = J pa 2 da. The global entropy maxi- 
mum of this problem is the distribution ([20]) . Then, we 
can express the entropy ([2]) in terms of uo and uo 2 by sub- 
stituting the optimal distribution ([20]) in Eq. (j2j). This 
gives the functional ([23]) . 

Second step: The maximization problem (|A1|) is now 
equivalent to 

max {5[SJ,cj2] | Vp] = T, Ep] = E, 

UJ,0U 2 

rf<V] = rf 9 }, (A2) 

where Sp, uo 2 ] is given by Eq. ([23]) . To solve that prob- 
lem, we now maximize Sp, uo 2 ) at fixed E, T, T^' 9 ' and 
a fixed vorticity profile uo(r). Since the specification of 
cJ(r) determines E, T and Y 2 9 ' , this is equivalent to max- 
imizing Sp, uo 2 ] at fixed vorticity profile cJ(r) and fixed 
J U02 dr = T2' 9 ' — Y 2 9 ' ■ The global entropy maximum of 
this problem is uu2(r) = ^2 = (^2 9 ~ F 2 g ')/A, where A 
is the domain area. Then, we can express the entropy 
([23]) in terms of uo by substituting the foregoing solution 
in Eq. ([23]) . This gives the functional 

S[ZJ] = | In n 2 = ^ In (r^- - r^- [ID]) , (A3) 

up to an additional constant. 

Third step: The maximization problem (|A2[) is now 
equivalent to 

max {Sp} I Tp] = T, Ep] = E}. (A4) 



where Sp] is given by Eq. (|A3|) . Since ln(x) is a mono- 
tonically increasing function, we can finally remark that 
this maximization problem is equivalent to 

min {T c 2 9 -p] I r[cJ]=T, Ep]=E}. (A5) 

In conclusion, the maximization of entropy at fixed circu- 
lation, energy and microscopic enstrophy ()A1|) and (jA2|) 
are equivalent to the minimization of coarse-grained en- 
strophy at fixed circulation and energy ([A5|) [32|. 

Remark: Writing the variational problem associated to 
(|A4[) in the form SS — f35 E — aST = (first variations), 
we recover Eq. ([27]) . 



Appendix B: Equivalence between the stability 
criteria ([38]) and ([58]) 

In Sec. IIIIBI we have shown the equivalence of ([36]) 
and ([52]) for global maximization. Another proof is given 
by Ellis et ai [HI by using large deviations technics (they 
show that the probability density of the coarse-grained 
vorticity Pp] is given by a Cramer formula involving 
the entropy ([44]) - ([47]) . so that the most probable coarse- 
grained vorticity solves the maximization problem ([52]) ). 
In this Appendix, we show the equivalence of ([36]) and 
([52]) for local maximization, i.e. p(r, a) is a (local) max- 
imum of S x [p] at fixed E, T and normalization if, and 
only if, the corresponding coarse-grained vorticity cJ(r) 
is a (local) maximum of Sp] at fixed E and T. To that 
purpose, using a method sketched in [52|, we show the 
equivalence between the stability criteria ([38]) and ([58]) . 
Another proof is given by Bouchet [22| (see also Appendix 
G of Chavanis |23|) by using an orthogonal decomposi- 
tion of the perturbation [17j. 

We shall determine the perturbation 5p*(r,a) that 
maximizes S 2 J[Sp] given by Eq. ([38]) with the constraints 
5uo = f Spa da and J Spda = 0, where Suo(r) is prescribed 
(assumed to conserve energy and circulation at first or- 
der). Since the specification of Suo determines Sip, hence 
the second integral in Eq. ([38]) . we can write the varia- 
tional problem in the form 

S ^-^ J ^^drda^j - J X(r)S Q Spada^j dr 

- J ((r)S (^J Spda^j dr = 0, (Bl) 

where A(r) and £(r) are Lagrange multipliers. This gives 
5p, = -p(r,(T)(A(r)a + C(r)), (B2) 

and it is a global maximum of S 2 J[Sp] with the previ- 
ous constraints since S 2 (S 2 J) = — J 6<y6 2 p J drda < (the 
constraints are linear in Sp so their second variations 
vanish). The Lagrange multipliers are determined from 
the constraints Suo = j Spa da and J Spda = yielding 
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Suj = — Xuo 2 — (uj and = — Xuj — (. Therefore, the optimal 
perturbation (|B2|) can finally be written 



£p* = —p(o- - uj). 
0J2 



(B3) 



Since it maximizes 5 2 J[5p] 1 we have S 2 J[Sp] < 5 2 J[5p*]. 
Explicating 5 2 J[5p*] using Eqs. ([38]) and (jB3[) . we obtain 

£ 2 J[£p] < - ^^dr- J SuS^dr. (B4) 

Finally, using Eq. ([5L]) . the foregoing inequality can be 
rewritten 



£ 2 J[£p] < ~j C" (uj)(5uj) 2 dv 
J Suj5t/jdr = 6 2 J[5uj}, 



(B5) 



where the r.h.s. is precisely the functional appearing in 
Eq. ([58]) . Furthermore, there is equality in Eq. ()B5|) 
iff £p = 5p*. This proves that the stability criteria ([38]) 
and ([58]) are equivalent. Indeed: (i) if inequality ([58]) is 
fulfilled for all perturbations 5uJ that conserve circulation 
and energy at first order, then according to Eq. (jB5|) , we 
know that inequality ([38]) is fulfilled for all perturbations 
5p that conserve circulation, energy, and normalization 
at first order; (ii) if there exists a perturbation SuJ c that 
makes S 2 J[Suj] > 0, then the perturbation 5p c given by 
Eq. ()B3|) with Sou = SuJ c makes S 2 J[Sp] > 0. In conclu- 
sion, the stability criteria ([38]) and ([58]) are equivalent. 



Combining these two relations, we obtain 



dp 

dt 



u • Vp = p(a- id)C"{id) I 



' duj 
\~dt 



(C5) 



Substituting Eq. (|6Zl in Eq. ([C5]), we find that 



dp „ 



/o(<t - w)C"(w)V • <^ D 



,(C6) 



Finally, using Eqs. (]C4|) and ([5T]h the foregoing equation 
can be rearranged in the form 



dp 
dt 



u - Vp 

C^2 



p((T - 07) 



Vp + /3(t)p(cr - cJ) 



(C7) 



Under this form, we see some analogies (but also crucial 
differences) with the relaxation equation ([T5]) related to 
the MRS approach. Finally, in the case of a Gaussian 
prior, cj2(r,t) = ^2 is constant and the foregoing equa- 
tion reduces to 



dp 
dt 



+ u-Vp 

= p(o — uj) V 



D 



p(a-uj) 



Vp + P(i)p(a -uJ)Vip 



Appendix C: The equation for the vorticity 
distribution 



In the parametrization of Sec. IIV Al the mean vorticity 
cJ(r, t) evolves according to Eqs. ([67]) -(|68 ]) and the vortic- 
ity distribution p(r, cr, t) is then given by Eqs. ([62]) - ([64]) . 
It can be of interest to determine the relaxation equation 
satisfied by p(r, cr, t). According to Eq. ([62]) . we have 



(Cl) 



Lap = -a^> + In x(cr) - lnx($), 



where $(r,£) is related to uJ(r,£) according to Eqs. 
([64]) . Taking the derivative with respect to time and 
using Eq. ([M]) . we obtain 



<91np _ N <9$ 

__ = _ ((7 _ w) _ 



Then, using Eq. ([46]) . we get 

Similarly, we have 

Vlnp= (a-w)C"(w)Vw. 



(C2) 



(C3) 



(C4) 



Appendix D: Derivation of the relaxation equation 



We write the dynamical equation as 

duj ^ 

— + u • VuJ = X, 

at 



(Dl) 



where X is an unknown quantity to be chosen so as to 
increase S[uj] while conserving E and Y. The time varia- 
tions of S are given by 

S = - J C\uJ)Xdr. (D2) 

On the other hand, the time variations of E and Y are 

E = I Xijjdr = 0, (D3) 

Xdr = 0. (D4) 



Following the Maximum Entropy Production Princi- 
ple, we maximize S with E = Y = and the additional 
constraint 



X 2 



<C(r,t). 



(D5) 
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The variational principle can be written in the form 
SS - 0(t)5E - a(t)St - J p^-j5 (^j dr = 0, 

(D6) 

and we obtain the following current 

X = —D (C\uj) + P(t)i/; + a(t)) . (D7) 

Substituting Eq. (JDZ)) in Eq. (jE]), we obtain Eq. ([75]) . 
The Lagrange multipliers f3(t) and a(t) evolve so as to 
satisfy the constraints. Substituting Eq. (|D7|) in Eqs. 



(|D3|) and (|D4[) . we obtain the algebraic equations ([76]) - 
([77]) . Finally, the H-theorem ([79]) can be derived as fol- 
lows. Substituting Eq. ([D7]) in Eq. (fD2l) . we get 

S= J X + /3(t)^ + a(*)^ dr. (D8) 

Then, using the constraints (|D3|) - (jD4|) we obtain 

/X 2 
— dr>0. (D9) 
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